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Ab stract — A statistical model is given for soil development relating meteoroid 
impacts on the moon to cosmic-ray-produced isotopes in the soil. By means of this 
model, the average lunar mass loss rate during the past 1.4 aeons is determined to 
be 170 g/cm aeon and the soil mixing rate to be ~200 cm/aeon from the Gd isotope 
data for the Apollo 15 and 16 drill stems. The isotope data also restrict the time 
variation of the meteoroid flux during the past 1.4 aeons. 


INTRODUCTION 

The lunar regolith is very complicated. Each gram of soil contains millions of 
separate grains, each with a distinct history. The histories of a few of the larger 
grains (rock chips) have been traced by isotope analysis; however, it would be 
impossible to trace the separate histories of all the grains that have been recovered. 
Some similarities should exist in the grain histories because the cosmic-ray bom- 
bardment is uniform over the moon, and to a good approximation, the bombardment 
rate by sufficiently small meteoroids can also be considered uniform over the moon 
during aeons. From these two uniformity conditions and from neutron isotope data in 
a statistical model of soil development described by an integral equation with rigid 
boundary conditions, we found that the moon has been losing mass during the past 
aeon at a rate of approximately 80 g/cm aeon (Fireman, 1974). We shall examine the 
consequences of these two uniformity conditions for a range of boundary conditions and 
also consider how the conditions change when the meteoroid flux increased in the past. 
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Cosmic -ray -produced isotopes have been measured as a function of depth to almost 
3 m deep in lunar soil. It is well known that the cosmic-ray exposure age is obtained 
from a stable and radioactive cosmic-ray-produced isotope when the relative produc- 
tion rates are known. All soil samples contain large amounts of solar wind. Since 
appreciable amounts of stable rare-gas isotopes are in solar wind, it is difficult to 
measure the spallation-produced stable rare-gas isotopes normally used for exposure- 
age estimates. In spite of this difficulty, Bogard et al . (1973) and Pepin et al . (1974) 
have measured the concentrations of spallation-produced stable isotopes as a function 
of depth. Their interesting result is that the depth variation of the spallation-produced 
stable rare-gas isotopes resembles that of neutron-produced isotopes rather than that 
of spallation radioactivities. 
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The spallation-produced radioactivities H and Ar (Fireman et al. , 1973) have 

37 235 

maxima near the surface, while the neutron-produced isotopes Ar' , U (n,f), and 
B*°(n,a) (Stoenner et al . , 1973; Stoenner and Davis, 1974; Woolum and Burnett, 

1974a, b) have maxima near 100-cm depth. The neutron-produced stable isotopes of 
Gd and Sm were measured and have maxima at depths greater than 100 cm (Russ et al ., 
1972; Russ, 1973). 

If no soil motion occurred during cosmic-ray irradiation, then the depth variation 
of the spallation-produced stable rare-gas isotopes would be identical to correspond- 
ingly produced radioactive isotopes, except for fluctuations near the surface caused 
by solar flares, and the depth variation of the neutron-produced isotopes of Gd and 
Sm could be calculated from cosmic-ray theory (Lingenfelter et al . , 1972). 

The depth variations of the spallation stable rare-gas isotopes and the neutron- 
produced isotopes of Gd and Sm are not as would be expected if there were no soil 
movement during the irradiation. Russ et al . (1972), Russ (1973), and Bogard et al . 
(1973) have proposed accretionary-type soil models that account for the differences 
in depth variation. Accretionary models necessitate that a large amount of irradiated 
soil be buried below 3 m, if the moon has been bombarded by cosmic rays for 4. 5 aeons. 
This required thickness is at least a factor of 2 larger than the soil -thickness estimates 
from crater data (Shoemaker et al . , 1970) and from active seismic investigations 
(Watkins and Kovacks, 1973), which range between 3 and 12 m. If the ancient soil has 


2 



escaped from the moon, then a regolith thickness even as small as 250 cm is con- 
sistent with a cosmic-ray irradiation for 4. 5 aeons (Fireman, 1974). Accretionaiy- 
type soil models are also inconsistent with the type of soil motion expected from 
meteoroid impacts. Since soil removed from one site and deposited at another by 
meteoroid impacts is generally deposited at a shallower depth, irradiated soil would 
be enhanced near the surface. We mention three other aspects of the statistical 
model: (1) A statistical model can be checked by isotope data at two different sites. 

The accretionary models consider each site to be so unique that there is little 
possibility for checking one of these models from isotope data at a different site. 

(2) With the statistical model, the average mass loss rate of the moon can be estimated. 

(3) Also, with the statistical model, isotope data can be related to meteor impact data. 

A GENERAL COSMIC-RAY REQUIREMENT FOR SOIL MODELS 

The stable cosmic- ray-produced isotopes at all six Apollo sites require that 
irradiated soil be diluted by unirradiated (less irradiated) soil during times of the 
order of aeons. Such dilution is far greater than the implantation of meteors with 
low exposure ages would have caused, since the chemical composition of the soil 
limits the amount of meteoric material to less than 2% (Laul et al . , 1971). 

Cosmic -ray -produced isotope data are most relevant for times between 0. 1 and 
4 aeons and depths between ~10 and 300 cm. On the other hand, solar- wind and 
solar-flare-track data are relevant for much smaller times and depths. All soil 
grains were exposed to solar wind for approximately a thousand years, and a signifi- 
cant fraction of grains were exposed to solar flares for a million years. Models for 
soil development from cosmic-ray data can be made consistent with these solar-wind 
and flare data by requiring the soil grains to spend a short time during their history 
near the surface. 

Also, occasionally, soil layers approximately 1 cm thick have a somewhat differ- 
ent color and texture than does the soil above and below. The most distinctive of 
these layers is at ~13-em depth for an Apollo 12 location (Laul et al . , 1971). This 
layer was essentially undisturbed for approximately 10 m.y. Models for soil develop- 
ment from cosmic-ray data should therefore permit a 1-cm layer of material at 
approximately 13-cm depth to persist undisturbed for more than 10 m.y. 
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From the point of view of cosmic rays, there is nothing unique about the six 
Apollo sites. If irradiated soil is diluted by unirradiated soil at these sites, it is 
natural to expect irradiated material to be diluted everywhere on the moon. If the 
cosmic-ray bombardment rate determined from lunar material is uniform over the 
moon, then the rate of dilution of irradiated material is also uniform. This require- 
ment is certainly true when soil dilution is caused by small meteor impacts during a 
long time. We impose this requirement on our soil model for times between 0. 1 and 
4 aeons and depths between 10 and 300 cm. We examine those soil models that satisfy 
this uniformity requirement for consistency with the Apollo 15 and 16 isotope data. 


The uniformity condition requires a relation between the dilution rate q(x) and 
the mass escape rate of soil from the moon, m • 

d 

m e = p / q(x) dx (g/cm 2 aeon) 

0 


1 

P 


dm 



(i) 


where q(x) is the fraction of mass at depth x per aeon replaced by unirradiated material, 
p is the density of the soil, and d is the depth of undisturbed soil. If all the irradiated 
grains that are replaced by unirradiated grains return to the moon, then the moon 
would have more irradiated material per unit area everywhere else than at the site 
under consideration. Relation (1) eliminates the uniqueness of the site. 

The dilution of irradiated soil affects the relation between the cosmic-ray isotope 
exposure age and the isotope production rate. The modification in the exposure age is 
discussed in the next section. 


THE SOIL MODEL 


Lunar soil has developed by meteoroid impacts. Theoretically, the meteoroid 
impact parameters can be calculated from the soil movement and mass change if they 
are known. The cosmic-ray-produced isotopes in the soil and, especially, their depth 
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variations also depend on the soil movement and its mass change. Since these isotope 
measurements are quite accurate, it is desirable to relate them to properties of the 
soil movement and the mass change for the moon. Our model formulates relations of 
this type. 

Meteoroid impacts below a certain size range can be regarded as uniform over 
the moon for a time scale of the order of aeons. Cosmic-ray-produced isotopes have 
been measured to depths of nearly 3 m and show approximately smooth variations over 
this depth range. An irregularity exists in the Gd isotope measurements of the 
Apollo 16 drill stem in the neighborhood of 200-cm depth (Russ, 1973), but it can be 
attributed to a difficulty that occurred during the sample collection. Even if this 
irregularity is real, statistical models could encompass it. 

Isotope measurements have been made at ^30-cm -depth intervals. We divide 
meteoroid impacts into two classes: those that remove or deposit less than ~30 cm of 
soil from a site, and those that remove or deposit more than m. The small impacts 
are assumed to be uniform over the moon during aeons; the large impacts are not. 
Intermediate impacts occasionally cause an irregularity in the isotope contents. 

Meteoroid impacts affect the soil by removing and replacing surface material 

and by mixing surface material to depth. Soil layers can persist for long times even 

when the soil is uniformly mixed to a depth of 3 m in an aeon. If the meteoroid flux 

were constant in time, soil would then be uniformly mixed to 3-cm depth in 10 m.y. , 

and a layer at 13-cm depth would be preserved for 40 m.y. Laul et al. (1971) have 
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observed a soil layer at 13-cm depth that has an age of 10 to 10 yr. If the meteoroid 
flux increased in the past, a soil layer at 13-cm depth could be preserved longer. 

A statistical soil model was mathematically formulated by an integral equation 
with stringent boundary conditions (Fireman, 1974). Here, we describe the same 
model in physical terms and relax the boundary conditions. 

If the cosmic-ray-produced isotope clock at a site s is reset by a major local 
impact at time T g , then the exposure ages T(x) as a function of depth for a constant 
meteoroid flux are 


5 



-T q(x) T f -T q(x)1 

T(x) = T g e + \l - e J for T g q(x) < 1 , (2a) 

-T q(x) , [ -T q(x)"| 

T(x)=T g e + 2q^T) L 1 " e -1 for T g q(x)^l , (2b) 

if q(x) is the dilution rate caused by minor impacts. The replacement time for soil at 
depth x is l/q(x). If the meteoroid flux increased in the past, then the exposure ages 
would be approximately 

-T q(x) T f -T q(x)l 

T(x) = T g e S + -g Ll - e 5 J for T g q(x) < 1 , (3a) 

-T q(xj if -T q(xjl 

T(x)=T g e + 2 ^ll-e J for T g q(x)>l , (3b) 


where q(x) is the time-averaged soil dilution rate during T . If the meteoroid flux 

2 R t .w 2. b T g 

increased exponentially in the past as e , then q(x) is [q^(x)/2. 6 T g j (e - 1), 

where q^(x) is the present dilution rate. Equations (2) and (3) give the relation between 

the exposure age and the soil replacement parameters, T and q(x). 

s 


There are also relations between the exposure age T(x) and isotope production 
rates R(x). If soil is not mixed below depth d, then 


T(x) _ — — 
( ) Rg(X) 


(Gd 158 /Gd 157 ) - (Gd 158 /Gd 157 ). 


[1 + (Gd 158 /Gd 157 )] R 


for 


x > d 


158 


(x) 


(4a) 
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He is a frequently measured rare-gas spallation isotope. 


~ .158 . „ 

Gd IQ Q WPkll — 


measured neutron-produced isotope in lunar soils, and an initial ratio (Gd 158 /Gd 157 ) i 
can be estimated. The relation between spallation rare-gas isotopes and exposure age 
is well known; that between the Gd isotope and exposure age was derived by Eugster 
et al. (1970). If soil is uniformly mixed to depth d, then 
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T(x) = 


for x < d . 


He 3 (Gd i58 /Gd 157 ) - (Gd 158 /Gd 157 ). 

d - d 

(1/d) J E 3 (x) [1 + <Gd l58 /Gd 157 )] (1/d) J R(x) dx 

0 0 


(4b) 


We assume that soil is uniformly mixed from the surface to depth d, where the 
dilution rate is zero, or to bedrock at a rate C, where 

d = CT . (5) 


If exposure-age relations (3) are substituted into equations (4), the soil replace- 
ment and mixing parameters are related to the cosmic-ray-produced isotopes. For 
example, equation (3a) is equated to (4b) to get 


T 

s 
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(Gd 158 /Gd 157 ) - (Gd l58 /Gd 157 ). 

- 

[1 + (Gd 158 /Gd 157 )] (1/d) j R 158 (x) dx 


when x < d and T q(x) £ 1; and equation (3b) is equated to (4b) to get 

S 


"T <i( x ) j 

T e s +x= 

s 2q(x) 


Q Q 

r -T q(x)1 Ar 

[l-e S J = 


(Gd 158 /Gd 157 ) - (Gd^/Gd 157 ^ 


J R 38 (x) dx [1 + <Gd 158 /Gd 157 )] (1/d) J R 15g (x) dx 


( 7 ) 


when x < d and T o q(x) s l. 


We calculate the soil parameters q(x) and C from Gd data at the Apollo 16 site 
with maximal and minimal values for T 10 , the time of the last major impact, and then 
examine how well Gd ratios calculated with the same q(x) and C but with a different 
time, Tjg, fit the Gd measurements (Russ et al . , 1972) for the Apollo 15 site. 
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The maximum T^ g is 4.0 aeons, the solidification age of the rocks at the Apollo 16 
site. Production rates R^ g g (x) have been calculated by both Lingenfelter et al. (1972) 
and Kornblum and Fireman (1974), and the two calculations agree to better than a factor 
of 2. The values used in Tables 1 and 2 are essentially those given by Lingenfelter et al. 
(1972). The dilution rates q(x) are calculated with either equation (6) or (7) from the 
Apollo 16 Gd data and are given in Table 1. Two cases are considered: bedrock at 
250-cm depth and a 750-cm depth for undisturbed soil obtained from setting 4 q(250) - 1. 

The minimum T ig is^l„40 aeons, determined from equation (6) by setting q(dj = 0 
and by maximizing (1/d) J R 1 _ e (x) dx. This minimal impact time, which occurs when 

Q -Lt)o i cq i(--7 

d is approximately 250 cm, depends on the initial Gd /Gd ratio. Reasonable 

limits for (Gd /Gd ). are the terrestrial value, 1. 587, and the lowest lunar value, 

1.592 (Russ et al . , 1972). With the initial ratio of 1.592, is 1.40 aeons. The 

corresponding dilution rates are given in Table 1. The bottom row in Table 1 gives 

the lunar mass escape rates corresponding to each of the calculated dilution rates. 

The Gd data at the Apollo 16 site impose three conditions on our model: (1) The time 

of the last major impact at the Apollo 16 site was between 1,40 and 4. 0 aeons ago. 

(2) Soil mixes from the surface down at a rate between 60 and 180 cm/aeon. (3) The 

o 

average mass escape rate of lunar soil from the moon is 171 g/cm aeon for 1.4 aeons 

2 

and ranges from 236 to 316 g/cm aeon for 4. 0 aeons. 

We next examine whether any soil mixing and mass escape rates obtained from 
the Apollo 16 Gd data are consistent with the Apollo 15 Gd data. 

SOIL MODEL APPLIED AT THE APOLLO 15 SITE 

For a constant meteoroid flux, the soil dilution and mixing rates at the Apollo 15 
site are the same as those at the Apollo 16 site because of the uniformity require- 
ment. For a meteoroid flux that increased exponentially in the past as e““ , which 

is the largest proposed increase (Fechtig, 1971), the uniformity requirement relates 
the time-averaged dilution and mixing rates at the sites by 

2 * 6T 15 2 * 6T 15 

q 15 (X) = ~ 27 & T , ~ q 16< x) 311(3 C 15 (X)=: 2. 6 T . °16 (x) * 

e 16 - 1 e 16 - 1 

( 8 ) 
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From equation (4a), the exposure age below the mixing depth d for the Apollo 15 
site is 


(Gd 158 /Gd 157 ) - (Gd i58 /Gd 157 ). 

15 { } [1 + (Gd 158 /Gd 157 )J [0.70R(x)l * 

158 

since the Gd production rate at that site is approximately 0.70 times that at 
Apollo 16 (Lingenfelter et al . , 1972). From Russ et al . ’s (1972) isotope ratios, the 
time of the last major impact at the Apollo 15 site (exposure age at 250-cm depth) is 

1.596 - (Gd 158 /Gd 15? ). 

T (250) = . (10) 

a 2.596 (0.28 X 10 4 ) 


Unless the initial Gd isotope ratio at the Apollo 15 site differs from that at Apollo 16, 
the time of the last impact at the Apollo 15 site ranges from 0. 54 to 1. 20 aeons. 

Certain values for the soil mixing and dilution rates consistent with the Apollo 16 
Gd data are inconsistent with the Apollo 15 Gd data. Our aim is to find a soil mixing 
rate and a mass escape rate consistent with the data at both sites or to show that none 
exists. We first assume a time-constant meteoroid flux and later assume a meteoroid 
flux that increased in the past. We require for consistency a difference of less than 

7 f-Q 7 r»T 

0. 001 between the calculated and the measured Gd /Gd ' ratios. 

1 1 

The Gd /Gd ratios, calculated with the dilution rates given in Table 1 for 
which (Gd 158 /Gd 157 ) i = 1. 592, are plotted as a function of depth in Fig. 1. When the 
meteoroid flux is constant, the dilution rate is identical for the two sites even though 
the major impact times at the two differ. Figure 1 shows that the soil mixing rate 

O 

of 250 cm/1.4 aeons, or 178 cm/aeon, and the mass escape rate of 171 g/cm aeon 
are consistent and that the other soil mixing and mass escape rates given in Table 1 
are not. 

The question arises as to why, with an identical statistical model, we obtain 171 

2 o 

g/cm aeon for the mass escape rate here and 80 g/cm 4 aeon in Firemen (1974). The 
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primary reason is that we use the neutron capture rates calculated by Lingenfelter 
etal. (1972) here and those calculated by Kornblum and Fireman (1974) there. This 
accounts for 70% of the difference. The second reason is that here we use a 
(Gd^/Gd^ 7 ). of 1. 592 for Apollo 16 soil in order to make Apollo 15 and 16 have the 
same initial ratio, while Fireman (1974) employed an initial ratio of 1. 5866 for Apollo 
16 soil and 1. 592 for Apollo 15. Finally, there is a slight difference in the soil mixing: 
Here, we have the soil mix from the surface with a rate C; Fireman (1974) required a 
square-wave type of soil mixing at all depths. 

It is well established from crater information that the flux of large meteoroids 

increased greatly between 2. 0 and 4. 0 aeons ago, but it is not certain whether there 

has been much change in the flux of small meteoroids during the past aeon. Fechtig 

(1971) has presented evidence that the meteoroid flux increased for all times in the 

past as e"* , where t is in aeons. It is of interest to examine how this meteoroid- 

flux change during the past 1.4 aeons would affect our soil model. With T,. = 1.4 aeons 

lo 

and T^ 5 = 0. 54 aeon, the time-averaged soil mixing rate at the Apollo 15 site is 

18 cm/aeon (approximately 10% of the Apollo 16 rate). Since the surface soil is 

mixed to a depth of only 9. 7 cm, the neutron capture rate is so low that the Gd^ 58 /Gd* 

ratio at x = 0 is not significantly altered from 1. 592. The soil impact parameters are 

not consistent with the Gd isotope data if the meteoroid flux increased for all times 
2 6 1 

in the past as e . Consistency is achieved, however, if the flux increase started 
1 aeon ago. 
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RATIOS IN APOLLO 15 SOIL 
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T |6 = 4.0 aeon, 250-cm DEPTH BEDROCK 



DEPTH (cm) 


Fig. 1. 


Comparison of calculated and measured Gd 


158 /Gd 15? 


ratios. 
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Table 1. Dilution rates q(x) 


158 157 

and mass escape rates determined from Apollo 16 soil Gd /Gd data (Russ, 1973). 





A/(l + (Gd 158 /Gd 157 )] 



A/[l + (Gd 158 /Gd l " j7 >l 


Depth 

(cm) 

R 158 W 

A 

rrt 158 

J -592 

Gd ilU 

/ 250 

[mf 

0 

R 158 WdX ) 

q(x) 1 

(iieon^j 

q(x) J 

, -J, 

<aeon ) 

/ 750 \ T 

\75o/ R 158 (X>dX ) 

0 

A >. 

q(x) 

(aeon" 1 ) 

0 

~0.1 X10" 2 

0.011 


0.69 

0.85 

1.05 

1.38 

0. 51 

35 

0* ** 60 X lo" 2 

0.014 


0.82 

0.77 

0.87 

1.64 

0.44 

70 

0. S7X 10” 2 

0.017 


0. 96 

0.67 1 

0.71 

1.92 

0. 38 

100 

1.00 X 10~ 2 

0.019 


1.07 

0.62 

0.45 

2. 14 

0. 33 

135 

0* 95 X 10' 2 

0.020 


1.13 

0. 60 

0.35 

2. 26 

0.32 

170 

0* 75 X Kf 2 

0.021 


1.18 

0.58 

0.27 

2.36 

0.30 

200 

0. 52 X 10" 2 

0. 022 


1.26 | 

0.54 

0. 17 

2. 52 

0.28 

250 

0. 38 X 10" 2 

0.25 


1.40 

0. 50 

0 

2.80 

0.25 

mass escape rate <m e ) 

250 

- P f Q(x) d* 

"0 

1< x i) 

2 

236 g/cm aeon 

1 

1 

9 

171 g/cm" aeon 

250 750 

p | q(xj dx 4- p j q(x) dx 
0 *250 

r 750 1 

129 4 P J* <l<x> dx g/cm 2 aeon 
L 250 


* 250 750 

250 / R 158 < x) * = °* 68 X 10 ‘ 2; / R 158 <*> ^ 34 X 10 ' 2 

0 0 

V 10 - 4.0 aeon; bedrock is at 250 cm, and q(x) is determined from equation <7). 

■^T 1( , = 1.4 aeon; q(250) = 0, which means soil is undisturbed at x = 250 cm: q(xj is determined from equation (7) for x > 100 cm anti from equation (6) for x < 100 cm. 

**T. ^ t= 4. o aeon; bedrock is at 750 cm and q(x) is determined from equation (7). 
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Table 2. 


Calculated Gd 158 /Gd 157 


ratios for Apollo 15 soil with dilution rates q(x) from Table 


* 

1 . 


Depth 

(cm) 

R 158 W 

q(x) t 

(Gd 158 /Gd 157 ) t 

II 

d-h 

(Gd 158 /Gd 15 V 

** 

q(x) 

(Gd 158 /Gd 157 ) 

0 

-0.07 X 10" 2 

0.85 

1. 5974 

1.05 

1.5948 

0. 51 

1. 5979 

35 

0.42 X 10“ 2 

0.77 

1.5975 

0. 87 

1. 5967 

0. 44 

1. 5980 

70 

0.61 X 10“ 2 

0.67 

1.5976 

0.71 

1. 5992 

0.38 

1. 5980 

100 

0.70 X 10" 2 

0.62 

1. 5977 

0.45 

1.6005 

0. 33 

1.5981 

135 

0.67 X 10~ 2 

0.60 

1.5977 

0.35 

1.6007 

0.32 

1. 5981 

170 

0. 527 X 10" 2 

0.58 

1.5978 

0.27 

1.599 

0.30 

1. 5982 

200 

0.365 X 10 -2 

0.54 

1. 5978 

0. 17 

1.597 

0. 28 

1. 5982 

250 

0. 276 X 10" 2 

0.50 

1. 5979 

0 

1.596 

0.25 

1. 5982 


^Ratios calculated with T^^ = 0. 54 aeon and Gd lo8 /Gd 157 = 1. 592. 

^Since q(x) T.^ < 1 and bedrock is at 250-cm depth, the ratios are calculated with equation (6) and 
250 

ikf R 158 «*=°-^5X10- 2 . 

0 

^ Since q(x) T < 1 and mixed depth C T . = 34 cm, the ratios are calculated with equation (4a) for x s 35 cm; the 

15 lo 34 

if -2 

ratio is calculated at x = 0 with equation (6) and j R 15 g (x) dx = 0. 25 X 10 . 

Since q(x) T ^ < 1 and soil is undisturbed below 260-cm depth in 0. 54 aeon, the ratios are calculated with 
equation (6). 



